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Abstract 


A unified viscous theory of 2-D thin airfoils and 3-D thin wings is developed with 
numerical examples. The viscous theory of the load distribution is unique and tends to the 
classical inviscid result with Kutta condition in the high Reynolds number limit. A new 
theory of 2-D section induced drag is introduced with specific applications to three cases of 
interest: 1) constant angle of attack, 2) parabolic camber, and 3) a flapped airfoil. The 
first case is also extended to a profiled leading edge foil. The well-known drag due to 
absence of leading edge suction is derived from the viscous theory. It is independent of 
Reynolds number for zero thickness and varies inversely with the square root of the 
Reynolds number based on the leading edge radius for profiled sections. The role of 
turbulence in the section induced drag problem is discussed. A theory of minimum section 
induced drag is derived and applied. For low Reynolds number the minimum drag load 
tends to the constant angle of attack solution and for high Reynolds number to an 
approximation of the parabolic camber solution. The parabolic camber section induced 
drag is about 4 percent greater than the ideal minimum at high Reynolds number. Two new 
concepts, the viscous induced drag angle and the viscous induced separation potential are 
introduced. The separation potential is calculated for three 2-D cases and for a 3-D 
rectangular wing. The potential is calculated with input from a standard doublet lattice 
wing code without recourse to any boundary layer calculations. Separation is indicated in 
regions where it is observed experimentally. The classical induced drag is recovered in the 
3-D high Reynolds number limit with an additional contribution that is Reynolds number 
dependent. The 3-D viscous theory of minimum induced drag yields an equation for the 
optimal spanwise and chordwise load distribution. The design of optimal wing tip 
planforms and camber distributions is possible with the viscous 3-D wing theory. 
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NOTATION 


A 

a -> ^exp 

B 

b 

c 

cd,Cd 


cl, Cl 

^turb 

E 

El(. z) 

F{x), F{x,y) 
fix), f(x,y) 
Hix) 

K,E 

K 0 

Kd 

Ke 

Kl 

K v 

tix), £ix,y) 

L 

Uy) 

cn 

P 

QO’Ql 

R 

Re 

Rn 


aspect ratio, b 2 /S 

profiled airfoil drag parameter, see (2.2.55-57) 
vorticity stream function, see (3.1.8) 
wing span 

2- D airfoil chord or mean aerodynamic chord 
drag-due-to-lift coefficient (2-D, 3-D) 
minimum section drag coefficient 

lift coefficients (2-D & 3-D) 

drag due to turbulent enstrophy, see (2.2.58) 

see (2.2.39) 

exponential integral, see (3.1.26) 

universal minimum drag load, see (2.2.32) or (3.1.31) 

camber shape 

Heaviside step function 

elliptic integrals of the first and second kind 

modified Bessel function 

section induced drag kernel function, see (2.1.35) and Figure 2-3 or 
(3.1.32) and Figure 3-2 

viscous part of the 3-D induced drag kernel, see (3.1.24) 

load kernel function, see (2.1.16) on Figure 2-2 or (3.1.20) and 

Figure 3-11. 

3- D kernel for separation potential, see (3.1.37) and (3.2.3). 
load distributions 

lift 

span load distribution 
nose radius of profiled airfoil 
pressure 
source strengths 

Euclidean distance between two points 
Reynolds number — ref chord 
Reynolds number — ref nose radius 
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s 

T 


M,v,w 

Uoo 

V 

( x,y,z ) or (4 ,ti,C) 
XoOO, x e (y ) 


planform area 
Trefftz plane 

Cartesian velocity component 
See stream velocity 
volume of all space 
Cartesian coordinates 
wing leading and trailing edges 


a,a(x), a(x,y) 
a d 
a v 


m 

A(p 

A<p e 

r 

Y 

Yx.Yy 

P 

a 

ca 


angle of attack, negative slope of the camber shape 
viscous induced drag angle, see (2.1.31) 
separation potential, see (2.1.48) or (3.2.4), also Figure 2-12 
and Figure 3-3. 
delta function 

doublet strength, see (3.1.39) or (3.3.19) 

wake doublet strength or section load (3.3.44) 

total circulation 

circulation/unit length 

3-D circulation components, see (3.2.6) 

dynamic and kinematic viscosity 

density 

Reynolds number — ref quarter chord 
vorticity 


y2 
I I 

m 

(< 7 ) 


Laplace operator 

denotes absolute value 

integral of F over chord or planform 

vector quantity q 
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I. 


Introduction 


The work reported herein is a direct outgrowth of two previous studies, Yates 
(1980), and Yates and Donaldson (1986). The earlier work on 2-D viscous thin airfoil 
theory (also see Yates 1982 for 3-D wing results) focused on the direct role of viscosity in 
the generation of lift. The direct effect of the primary viscous boundary layer was taken out 
of the lift problem by considering the surface to move tangentially at the speed calculated 
with potential theory in its non-lifting attitude — the so-called "belt sander" model. We use 
this concept throughout the present study. When the airfoil is placed at a small angle of 
attack, the perturbation viscous boundary conditions induce an interactive vorticity field due 
to the surface pressure gradients. This field of local vorticity is responsible for the lift via 
the integral relation 


L 



S 


( 1 . 1 ) 


where S is the 2-D region of flow outside the foil. The complete theory of the 2-D foil 
including finite thickness and unsteady motion was documented in Yates (1980). A very 
important feature of the viscous wing theory is that the load distribution is uniquely 
determined without recourse to the empirical Kutta condition. Furthermore, the lift is 
established by the action of viscosity over the entire wetted surface and not merely the local 
region near the trailing edge. In this early study we calculated the lift as a function of 
Reynolds number. For Re «> , the lift tends to the Kutta Joukowski value 

c L = 27ta (1.2) 

in accordance with earlier results of Shen and Crimi (1965) with corrections of order 
(h/Re). Further calculations for finite trailing edge thickness indicate a proportional 
reduction of the lift. These and many other results were documented in our earlier work. 

The effort summarized in Yates and Donaldson (1986) was concerned primarily 
with the problem of drag-due-to-lift. An attempt was made to rationally assess various 
proposed and functioning devices for reducing this important drag component. To provide 
a consistent framework for this task, we revisited the theoretical foundations of the theory 
of drag. Many known representations for the drag were rederived using the momentum 
integral approach. By definition the drag can be expressed as the resultant of pressure and 



viscous forces on the surface or as a momentum flux integral in the Trefftz plane — all well- 
known results. The difficulty with using either of these formulas to attack the drag 
assessment problem is that two many flow variables are involved (velocity, pressure 
viscous stress, etc.). At the insistence of Donaldson, an all out effort was made to express 
the drag in terms of a single flow property. The trick for achieving this goal was to 
introduce the energy equation into the Trefftz plane momentum representation of the drag 
and then allow the Trefftz plane to recede to downstream infinity. This leads to a 
relationship for the power, U x times the Drag, as the total viscous dissipation in the fluid 
as required by thermodynamics. The crucial step was to show that in an incompressible 
fluid the total dissipation can be rewritten as the global enstrophy or that 

D = | \a)f cN (1.3) 

U oo rr 


where |i is the viscosity and a> is the vorticity. This is an exact relation that contains all 
components of drag by whatever engineering name we may give them. Since the lift can 
also be expressed entirely in terms of vorticity, 

L = pU ao \dS(D x y, (1.4) 

T 


where co* is the streamwise component of vorticity and T is the Trefftz plane, we have 
all components of the resultant force expressed in terms of vorticity alone. Many of the 
results in Yates and Donaldson (1986) are a direct result of these important formulas. The 
experimentally verified universal nature of the drag polar (e.g., see Rooney 1990) follows 
at once from the respective linear and quadratic dependence of the lift and drag on vorticity. 
Also, to the Authors knowledge, the dissipation or enstrophy formula for the drag is the 
only form that shows unequivocally that the drag is positive or in the direction of the 
external steady flow. 

During the course of the 1986 effort, the Author wrote a technical memorandum 
that outlined the development of a unified viscous approach to the theory of lift and drag 
(Yates 1986). The above formulas (1.3) and (1.4) together with the earlier work on 
viscous airfoil theory formed the core of the idea. We had seen that lift could be expressed 
entirely in terms of the interactive vorticity without reference to the vorticity in the primary 
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boundary layer - removed with the " belt sander" model. Even without the artiface of a 
moving surface we observed that the vorticity on a symmetric airfoil section at zero lift has 
odd parity (top to bottom) and thus integrates to zero under the spatial averaging operation 
(1.1). The interactive vorticity associated with the lift has even parity and yields a finite 
number upon integration. If we denote these two components symbolically by 
6> 0 and <w, we can write in 3-D flow 


(b = COq + 

|d)| 2 =|d)o| 2 + 2d)o-d) 1 +|d) 1 | 2 (1.5) 


But the product of an even and odd function is odd and thus integrates to zero. Thus 

L = pU„ j ds (o\y 

T 

and 

D = ^JN 2 ^ + 7t/N 2 ^ (1.6) 

u °° y U oo y 

The drag-due-to-lift is given by the second integral in (1.6). The idea that evolved was to 
use the interactive vorticity that we computed with the linear viscous theory in the quadratic 
expression for the drag, and thereby calculate the drag-due-to-lift from first principles. 
That is the primary subject of the present report. 

In Section II, we develop the basic viscous theory of the thin 2-D airfoil and present 
a collection of asymptotic and numerical results. The calculation of section drag-due-to-lift 
and the minimum section drag load distribution is a significant new contribution to the 
subject A detailed derivation of the section drag of a flat plate at angle of attack is given-a 
Reynolds number independent result that agrees with experiment. Also, a composite 
formula for the Reynolds dependent drag of a profiled section is suggested and compared 
with experimental data. Also, we introduce the concept of a viscous separation potential in 
Section II that follows from an expanded interpretation of our earlier viscous theory of lift. 
Preliminary examples indicate that this function could become a valuable design tool. 

In Section in, we develop the viscous theory for the 3-D wing with asymptotic and 
numerical results. The Reynolds dependent formula for drag-due-to-lift contains the 
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classical induced drag as well as the section drag-due-to-chordwise load variations of the 
2-D theory. The corresponding theory of minimum drag yields an integral equation for the 
optimal chordwise and spanwise load distribution. The theory of the 3-D separation 
potential for high Reynolds number deserves further development and application. A 
preliminary example for the rectangular wing indicates a large potential for separation at the 
wing tips as we might expect Also, the theory indicates that more pointed wing tips have a 
considerably reduced potential for separation. 

In conclusion, we remark that many of the results in this report were first outlined 
at an AIAA sponsored workshop on drag in which the Author participated as an invited 
lecturer (Yates 1990), under NASA Langley sponsorship. The AIAA notes should be used 
in conjunction with the present document. The author found that considerable scientific 
interest has been generated in the formula for drag (1.3). J. van der Vooren and J.W. 
Slooff presented work on the compressible counterpart of our 1986 formula while J.D. 
McLean presented detailed calculations of viscous drag that used the formula directly with 
input provided by boundary layer theory (see AIAA Professional Studies Series 1990). 
Earlier work of Greene (1988) on viscous induced drag was also discussed by Yates in the 
AIAA notes. A recent comprehensive discussion of many aspects of the compressible 
viscous theory may be found in Home et al. (1990). This important work was provided to 
the author by Home just prior to the AIAA drag workshop but was not included in the 
reference list. 
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II. Complete Viscous Theory of the 2-D Thin Airfoil 


2.1 Basic viscous theory 

The complete viscous theory of the two-dimensional zero-thickness airfoil is 
developed in detail. The theory is complete in the sense of Birkhoff s reversibility paradox: 

Any reversible (inviscid) theory of fluid dynamics is either incomplete, over 
determined or grossly incorrect, so far as its prediction of steady state lift 
and drag are concerned. 

Birkhoff (1955) 

By including viscosity at the outset we obtain a well posed problem for lift and drag that 
does not require auxiliary empirical conditions as does the inviscid theory. The theory of 
lift follows closely the earlier work of Yates (1980). The complete and powerful theory of 
section induced drag grew out of a study of Yates and Donaldson (1986). Part of this 
work was included in a course on drag and lift for the AIAA (see Yates 1990). The theory 
of the separation potential is completely new to the author's knowledge. For completeness, 
we develop the 2-D theory from scratch herein even though the theory of lift was first 
derived over a decade ago. Many of the ideas of the lift theory can be brought into much 
sharper focus when considered along with the drag problem. 

Consider a 2-D zero thickness airfoil of chord c fixed in a steady viscous 
incompressible flow ( U „ , p) along the x-axis. We scale all lengths by the semi-chord 
c/2, all velocities by (/» and pressure by P^» . The natural Reynolds number of the 
problem is referred to the quarter chord: 


<7 = 


UooC 

4v 


( 2 . 1 . 1 ) 


where v is the kinematic viscosity. The problem is illustrated in Figure 2-1. 
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Figure 2-1. The 2-D airfoil problem. 

The present theory is concerned solely with the prediction of lift and drag due to small 
perturbations y =f{x) of the surface C : y = 0 ±; W < 1. Following Yates (1980) we use 
a "belt sander " model to remove the mean boundary layer on C. Suppose that the upper 
and lower surfaces of C move at the free stream velocity. Then the exact viscous solution 
of the flow past C is the trivial potential flow solution. When we perturb C to the state 
y=f(x), the viscous problem for the perturbation velocity, pressure and vorticity 

becomes. 


du dv n 
dx dy 
dp du _ 
dx dx 


dp d\ 
dy dx 


1 dtp 
2 a dy 
1 da 
2a dx 


03 = 


dv 

Ik 


du 

% 


( 2 . 1 . 2 ) 
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with the viscous boundary conditions 


u = 0 and v = fix) = -a(x) on y = 0± 


(2.1.3) 


The homogeneous no-slip condition on the perturbation velocity follows from the use of 
the belt sander model and distinguishes the present theory from the conventional Oseen 
model (Shen and Crimi 1965, also Southwell and Squire 1933). By cross differentiating 
the momentum equations in (2.1.2) we obtain a useful alternative formulation in terms of 
pressure and vorticity: 


V 2 p = 0 


dp) 

dx 


— V 2 ft) 
2 a 


with the boundary conditions 


dp 

dx 


rw+f= 

dy 


1 d(0 
2a dy 

1 do) 

2 a dx 


on y = 0± 


(2.1.4) 


(2.1.5) 


These are the physical boundary conditions that relate the production of the linear 
interactive vorticity to gradients of the pressure that are driven by the geometric disturbance 
f(x). The vorticity formulation of the problem is fundamental since the lift and drag due to 
the camber disturbance f(x) can both be expressed in terms of to; i.e., 

Cl = KpUTc 

c ^WJTc = ht ldxdy (21 - 6) 



The above relations are exact as are the linear boundary conditions (2.1.5). The lift is 
proportional to the total vorticity or circulation that the foil can store — hence our 
interpretation that an airfoil is a vorticity capacitor (Yates 1990). The drag is the global 
integral of the square of the vorticity divided by the Reynolds number. For the present 
linear problem this exact result yields the section drag-due-to-lift. The profile drag has 
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been removed via the belt sander model. Our claim to completeness of the present theory 
follows directly from this elegant positive definite representation of the drag in terms of co. 

Integral representations 

The viscous perturbation problem can be recast as a problem for the single load 
function 


l(x) = p(x,0~)- p(x,Q + ) 


(2.1.7) 


The procedure is to use integral representations for p and co that follow from the 
fundamental solutions of (2.1.4). We have 


p = -^-jd% q 0 ($)InR 


with 


3 l 

m = ( 2 . 1 . 8 ) 


R = l(.x-$) 2 + y 2 P (2.1.9) 


where qo,qi are unknown distributions and Kq is the modified Bessel function 
(Abramowitz and Stegun 1964). The representation of the velocity components (u,v) 
follows from the momentum equations in (2.1.2) with the assumption that u,v vanish far 
upstream. Thus 


u = 


-P- 


iu^ 

2cr ' dy 

— oo 


V = 



0) 

2a 


( 2 . 1 . 10 ) 
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or with (2.1.8) 


u = -j-jdt q 0 ($)ln R--Lj-$dS q^^K^oR) 

v= j^^d^q^lnR + ^jj- q^)e a{x ~^K Q {oR) (2.1.11) 

Now compute the limit of u,v for y-»Q± and use the boundary conditions (2.1.3) 


_ 71 

lim u = +nqJx) ± — q. (*) = 0 
>-o± 2cr 

=*q l (x) = 2oq 0 (x) 


( 2 . 1 . 12 ) 


Also, 


l 

*J ± v = - \d$ q Q ^){tn\x- $\ + e a{x -^K 0 (<^x- = f'(x) (2.1.13) 

Finally, take the limit of p as y — > 0± and use (2.1.7) to get 

lim p = ±7tq 0 (x) => q 0 (x) = - ^ (2.1.14) 

>-°± 2 n 

Thus, we obtain an integral equation for the load function, 

1 i 

— J d$ t(S)K L (x - 4) = a(x) = -f(x) 

with 

^U ) = £[/nW + e m ^°(4<l)] 

The last result may be found in Yates (1980). The kernel function Ki of the load equation 
is an essential element of the viscous theory. It is plotted in Figure 2-2 versus the 
Reynolds scale X = ax and compared with the Cauchy singular inviscid result. No matter 


(2.1.15) 

(2.1.16) 
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how large the Reynolds number o, the logarithmic singularity of Ki is retained along 
with a more persistent downstream influence than upstream. This behavior is the reason 
for uniqueness of the viscous loads problem with a tendency to non-singular trailing edge 
behavior of the load in the high a limit. 



X - Reynolds scale 


Figure 2-2. Comparison of the viscous and inviscid 
load kernel functions, K^. 


With the solution of (2.1.15), the velocity, pressure and vorticity are given by the 
following integrals: 
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i a 1 

u = *(£)[*« * + e^-^^oCaR)] 

v = _ 2^^ Jd£4l)[*n*+e a(x '%o(<rt?)] 

p= "^i_J^ <(?)/ ' ,R 

. i 

® = — | ^ \d^l^)e o{x ~^K 0 {aR) (2.1.17) 


Section lift 

Next we use the above formula for to in the integral expression for the lift (2.1.6). 
First we derive an intermediate result for the circulation per unit length. Define 


oo 

Y(x) = - jdyo)(x t y) 

— oo 

= f it I d ^ ( t )e<T(X ~ 5) J^^o(oK) 

— 1 — oo 


(2.1.18) 


or using the known integral (Erdelyi 1954) 


we obtain 


]dyK 0 ^a(x 2 +y 2 f 2 y^e^ 

Y(x) = ■£■ t^)e a{x ~^ x ~^ 


(2.1.19) 


( 2 . 1 . 20 ) 


The section lift coefficient is given by 
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( 2 . 1 . 21 ) 


oo OO 

c L~~j C0 duty - J dx Y = T 

— oo — oo 

1 

= Jdx £(x) 

-l 


We further observe with (2.1.20) that 


y(x) = 2oe 2ax | d$ £(£)e 2a * 
-l 


= 0 


jc < — 1 
jc > 1 


( 2 . 1 . 22 ) 


The circulation decays exponentially on the scale of o upstream of the leading edge and is 
exactly zero downstream of the trailing edge. All of the vorticity that makes up the 
circulation y or the total circulation T is confined to the immediate neighborhood of the 
foil. This remarkable result is equally valid at high lift as is the general formula for lift 
(2.1.6). It is the viscous counterpart of the Kutta-Joukowski theorem. The ability of any 
foil configuration to develop lift is one and the same as its ability to store vorticity in the 
near field — hence our interpretation (Yates 1990) of an airfoil as a vorticity capacitor. We 
hasten to point out however, that the stored vorticity is continuously replenished via the 
steady state transport process in (2.1.4) and the production process in (2.1.5). At high 
angle of attack, these processes will in general become unstable with intermittent vortex 

shedding. 

Section induced drag 

The viscous formula (2.1.6) for the section induced drag is a very important new 
result that enables us to complete the theory of the viscous thin airfoil. With (2.1.6) and 
the integral representation of the vorticity (2.1.17) we can derive a useful practical formula 
for the section induced drag in terms of the load £(x). We substitute (2.1.17) into (2.1.6) 
and carry out the integration over the infinite domain. Thus 

l l 

c D =-?j \dt£(Z)\dil£(ri) x (2.1.23) 

2jt -l -1 

with 
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x- 


dxdy e^-^-»K 0 <M ( )K 0 (M n ) 


,<T(x-^+x-v), 


(2.1.24) 


where the subscripts (^,ri) on R = (x 2 +y 2 ) 1 ^ means that x is to be replaced by x-t, 
or *-r| respectively. The integration in (2.1.24) is completed as follows. 

First consider the integral 

t? = J_ oo ^y K 0 [a(a 2 + y 2 )% ]K 0 [a(b 2 + y 2 )% ] (2.1.25) 


Replace the second Bessel argument by its Fourier cosine representation (Erdelyi 1954) 


KQ[o(b 2 +y 2 ) 1/2 ] = jdt cos yt 


exp[-|6|(r 2 + o 2 )%] 
(t 2 + c r 2 )^ 


(2.1.26) 


then derive the Parseval relation for i}; 


* = J ^ cos ^ 0(A(a 2 + ^, 


o (r + « 


= - jK 0 (t)dt 
° <t(|<i|+|6|) 


(2.1.27) 


Thus, the function x can be written as 


d 2 




(2.1.28) 


<t(|j:-||+|x-tj|) 


or integrate by parts with respect to x to get 

n d 2 °° 

X= 2ad&j + (2.1.29) 


13 



Now we observe that is symmetric with respect to an interchange of Thus 

we assume without loss of generality that ^ > T| for the moment Then 




a 3&i] 

2 00 


n d 
+— 


a d%dr\ 


| dx (o{2x -$-»?)) 


k 


= 4r^7T- Jdf AT 0 (Osinhr 
= ^ ^[ k 0 ( o \$ - T7|)sinh <x(£ - r?)] 


(2.1.30) 


which is valid for all % and r| . Finally, we use (2.1.23) with the change of notation 
(^■q) — > (*,£) to obtain the final expression for the section induced drag 

l 

c D = jdxi(x)a d (x) 

-l 

a d (x) = ^ldSt{&^[K 0 ((fy-Z\)smho(x-Zj\ (2.1.31) 


where we interpret a d (x) as the viscous induced drag angle. It plays a role in the 2-D 
theory much like the induced angle of attack that is familiar in 3-D lifting line theory. 
Numerical examples of a d are presented in Section 2.3. The section drag is obviously 
quadratic in the section lift and is a function of Reynolds number a. There is no 
counterpart of this important result in the familiar inviscid formulation of the 2-D airfoil 
problem where the section drag is exactly zero. We make extensive use of this formula in 
the subsequent analysis. 

Minimum section induced drag 

With an explicit positive definite expression for the section drag-due-to-lift we 
can pose a variational problem for the minimum drag — a problem that is quite familiar in 
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3-D wing theory but heretofore non-existent in the 2-D theory. We seek a minimum of 
(2.1.31) with the constraint (2.1.21) on the total lift. Consider the functional 


l 

J[l] = D[(]-X[\dxl-c L ] 
-1 


(2.1.32) 


where X is a Lagrange multiplier. Now require J to be a minimum with respect to i(x) 
and X, 


5J = 0 and —— = 0 

dX 

The variation of J yields the Euler equation for fi. , 

-jdSt(0K„(x-4)=x 

71 ' 


(2.1.33) 


(2.1.34) 


with 


K °^~~dx sinh °*) 


The constraint on 12. is. 


| £(x)dx = 


-i 


C L 


(2.1.35) 


(2.1.36) 


and the minimum drag is obtained with (2.1.31) and (2.1.34), 


c d — 


CjX 

2 


(2.1.37) 


Since it is redundant to retain both X and cl, we define the universal load function F(x) 
with 
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(2.1.38) 



Then F(x ) is the solution of 


and 


or 


F($)K D {x-t) = 1 
l 

X = 2c L / \F(x)dx 

-1 


1 

cd = c l/ J ^00* 
-1 


^% = 1 / [F(x)dx 


(2.1.39) 


(2.1.40) 


(2.1.41) 


(2.1.42) 


*Cx) = cl • ^(*) / 


F(x)dr 


-l 


(2.1.43) 


An immediate result is evident if we compare the equation (2.1.39) for the 
minimum drag load with the induced drag angle in (2.1.31). The condition for minimum 
section drag is that the viscous induced drag angle be constant over the chord. Recall that 
this is also the condition for a minimum of the classical 3-D wing induced drag. The 
uniqueness of solution of the minimum drag equation (2.1.39) follows from the 
logarithmic behavior of the drag kernel K D near the origin. This kernel is plotted versus 
the Reynolds scale X = ox in Figure 2-3. 
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Figure 2-3. Section induced drag kernel, Kp. 


The camber shape and local angle of attack that produce the minimum drag load 
distribution are obtained with (2.1.15), 


/CO = - {|+««-njr 0 (o|*- $D] (2.1.44) 


with 


and 


G(x) = F(x) / \F(x)dx 
-1 


, , df 
a(x) = — j- 

dx 


(2.1.45) 


(2.1.46) 


The theory of minimum induced section drag has become a theory of the optimal load 
distribution. Our original viscous theory of the load distribution has become a theory of the 
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optimal camber shape. Recall that the classical theory of minimum induced drag of 3-D 
wings does not provide any information about the chordwise load distribution or shape. 
The above analysis is an important step in the direction of a true 3-D wing design code. 
Detailed calculations are presented in Section 2.3. 

The separation potential 

Since we have removed the primary boundary layer from the formulation of the 
viscous perturbation problem, it may seem pointless to discuss separation. However, the 
theory of the load function in (2.1.15) and (2.1.16) suggests an interesting approach to the 
problem. We write (2.1.15) in the form 

= «U) + «v(*) 

2n ll x ~£ 

with 

« V C 0 =- £jJ4 

The function a v 00 is a viscous load induced angle of attack or when multiplied by the 
free stream velocity it is a viscous load induced blowing velocity. We suppose that the real 
boundary layer that we removed must respond to this velocity. We show below with 
several examples that a v increases from leading edge to trailing edge where it becomes 
singular. For the class of load functions < 0 W that solve the corresponding inviscid 
problem 



(2.1.47) 


(2.1.48) 


the least singular a v (x) is obtained when 1 0 satisfies the Kutta condition at the trailing 
edge. This suggests that oc v Qt) may be used as an indicator of the potential for separation. 
We explore this idea in more detail in subsequent sections. 
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2.2 Asymptotic analysis 


In this section we examine the viscous theory more closely in the extreme limits of 
small and large quarter chord Reynolds number a. The analysis is of interest for several 
reasons. First we will show that the theory of the load function and the minimum drag load 
function is unique — an essential feature of any " complete " theory. Second, we gain further 
insight into the meaning of the separation potential. 

Low Reynolds number a « 1 

Consider the solution of the load equation (2.1.15) with the kernel function 
(2.1.16) when a tends to zero. For Lxl bounded, we first expand the kernel. For 
olxl « 1, we have 



+ e“^(o |*|)] 


+ In — + y + 1 


+ O(o 2 x in ojx|) 


( 2 . 2 . 1 ) 


where y is Euler's constant (= 0.5772). The integral equation for the load becomes 


^i|*-£| + iny + y+lj = a(;t) 


( 2 . 2 . 2 ) 


Within the class of load functions that admit a square root singularity at the leading and 
trailing edge, the last equation has a unique solution. This is a general feature of integral 
equations of the type (2.2.2) that have a log singular kernel function (e.g. see, Carrier, 
Krook and Pearson 1966). Since we know (see Fig. 2-2) that the logarithmic singularity is 
retained at all Reynolds numbers, we conclude that the solution of the load equation is 
unique for all a. For the same reason, the solution of the minimum drag loads problem is 
unique. 

We now demonstrate the solution of (2.2.2) for the special case of constant angle of 
attack. For then we can differentiate (2.2.2) to get 
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(2.2.3) 


q } m 


which has the eigensolution 


t(x) = 


Vi-? 


(2.2.4) 


Then substitute l(x) into (2.2.2) with oc(x) — a to get 


C = 
t(x) = 


2a 


a{ln 4 Jo - y - 1) 

2a 1 

a(tn 4/(7 — y - 1) Vl-x 2 


and 


c, = 


2tt a 


L c(tn4/G-Y-l) 


(2.2.5) 


( 2 . 2 . 6 ) 


The slope of the lift curve in low Reynolds number Stokes flow becomes infinite; i.e. 


dc L 


2n 


= 0 \ 


da (j{ln 4/<j — y — 1) ya In 4/c 


(2.2.7) 


We turn now to the minimum drag problem at low Reynolds number. First expand 
the drag kernel function Kp in (2.1.35) to get 

K D (x) = ^(K 0 (°W) sinh ax) 

= -<^ln\x\ + ln^ + y + \^ (2.2.9) 

which is identical to the low Reynolds number approximation of (2.2.1). Thus the 
universal minimum drag load function becomes 
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2 


1 


and 


F(x) = 


£(x) = 


(j(ln 4/<y- y-1) Vl-x 2 

c, 


Wl-JC 2 


— §■ = i / 



( 2 . 2 . 10 ) 

( 2 . 2 . 11 ) 

( 2 . 2 . 12 ) 


As we noted above, the minimum section drag problem has a unique solution since the 
kernel function Kq[x) has a logarithmic singularity. The induced drag number dc D / dcj_ 
tends to zero in the low Reynolds number limit. We will show presently that this number 
is a maximum for as 1. The minimum drag camber shape is the constant angle of attack 
fiiven by (2.2. 12), which tends to zero as o — > 0 for fixed cj_,. This is an example of 
two theorems that Batchelor (1967) attributes to Helmholtz:* 

1 . There cannot be more than one solution for the velocity distribution 
for flow in a given region with negligible inertia forces and consistent 
with prescribed values of the velocity vector at the boundary of the 
region, (including a hypothetical boundary at infinity when the fluid is 
of infinite extent). 

2. Flow with negligible inertia forces has a smaller total dissipation than 
any other incompressible flow on the same region with the same 
values of the velocity vector everywhere on the boundary of the 

region. 


Helmholtz 1868 
(see Batchelor 1967) 


*This author has tried unsuccessfully to obtain an extension of the second principle to higher Reynolds 
number via an enstrophy principle, (see Yates 1983). 
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High Reynolds number, o » 1 

The motivation and objective of the viscous theory is to obtain useful and practical 
results at high Reynolds number. By using the belt sander model we have formally 
removed any real Reynolds number associated with the primary boundary layer from the 
problem, (e.g., 8* or 0), even though a is one fourth of the usual chord Reynolds 
number. In the language of triple deck theory, (e.g. see Melnik (1978) or Brown & 
Stewartson (1970)), we have eliminated the middle deck in order to get at the inherent 
physics of the viscous origin of lift and section induced drag. Refinements that include 
Reynolds number effects of the middle deck and even geometric thickness can be 

considered later. 

First we write (2.1.15) in the form 

i-L U <(|)[<n|x-41+e^-«>jr 0 (<^-^l)] = «w < 2 ' 2 ' 13 > 

dx2n^ L 

and note that the contribution of the region lx - ^1 < X/c to the integral is small when 

1 « X « a. Thus, we can replace the Bessel function with its asymptotic expansion; 

i.e., 


e *°^~pojx( 


J-5- 

x> X/o 

V 2 ox 

= 0 

x<X/a 


(2.2.14) 


and so for e = X/c — > 0, we get 


-Lta-SS 

In 1 * r- 


m , i 


2jc ' x — £ 2-sJlna dx 


= a for a » 1 




(2.2.15) 
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where the first integral is a Cauchy principle value integral. The second term is the 
derivative of the Abel transform of the load. From the formulas (2.1.47) and (2.1.48) we 
recognize the second term as the viscous induced angle of attack; i.e.. 


and 


a v W = - 


1 d 
I'JlnO dx 




-1 


KS) 

Jx-1 



m 

x-4 


a(x) + a v (x) 


(2.2.16) 


(2.2.17) 


Constant alpha separation potential 

Our first use of oc v (a:) is to derive the asymptotic solution of (2.2.17) regarding 
a v as an error term for the corresponding inviscid problem 

-LjdS^ = a (2.2.18) 

27r _i X ~Z 


where for this example we assume a 

t(x) 

where cl is arbitrary. But with the viscous equation we can compute the error due to 
/( x). Thus, 


= constant. The solution of (2.2.18) is 


K 


c, - Inca 

. -Jl-x 2 , 


(2.2.19) 


l s r rff 


2-^lno dx ^ ^x-% 


i d ] 4 

f \ 

c L - 2 nal; 

Inline dx — % 

V *-« 2 J 


( 2 . 2 . 20 ) 


This error angle can be evaluated explicitly (see Appendix A). We have 
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with 



We recognize (2.2.23) as the solution of the inviscid problem with the auxiliary empirical 
statement of the Kutta condition. In a formal sense, we have derived the Kutta condition 
via the high Reynolds number viscous theory, a result first obtained by Shen and Crimi 
(1965) within the context of Oseen theory. 

Physically, the last result is much deeper than a simple empirical statement about the 
nature of the trailing edge flow. In Figure 2-4, we plot the normalized separation potential 

Vcr a v / cr for various values of C — . The reduced trailing edge singularity for 

2 na 

C L = 1 is clearly evident. It is important to note the monotonic growth from the 
2jccc 

leading edge to trailing edge where it is a maximum. This is our first example of the 
separation potential. The response of the upper surface primary boundary layer to this 
viscous blowing excitation will lead to separation at some distance from the trailing edge 
and this distance will increase with increased angle of attack. With repeated application of 


24 


the linear viscous theory and boundary layer analysis we could in principle develop a 
rational theory of high lift and stall. 



-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1 


x 

Figure 2-4. Separation potential and uniqueness of the load problem. 


Parabolic camber separation potential 

Another example of the separation potential is easily derived for the elliptic load 
distribution 

l(x) = ^^\- x 2 (2.2.25) 

K 

and corresponding parabolic camber distribution 
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(2.2.26) 


n n \ 


S 

X-$ 


_£i . 


;r 




The separation potential is given by 

c L d f 
av n^jlno dx 

which again can be evaluated in terms of elliptic integrals (Appendix A), 
parabolic camber 




(2.2.27) 
Thus for 


a v (x) = (K - 2£) (2.2.28) 

n^lna 

which we may compare with (2.2.24) for constant angle of attack. For x -> 1, we have 

a w (jc) = K a ~ constant 

and 

cc v (x) = — ^=K parabolic camber (2.2.29) 

n^no 


The trailing edge singularity in the parabolic camber potential is of the same logarithmic 
order but 4 times greater. Near the leading edge, x — » -1, we have 

aJx) = c L l\6^fno a = constant 

(2.2.30) 

= - cl / 2V net parabolic camber 

We interpret the negative value of a v as an indication of separation on the compression 
side of the leading edge. For the same lift coefficient and Reynolds number, the potential 
for separation of the parabolic cambered foil is 3 to 4 times greater than that of the foil at 
constant angle of attack. In Section 2.3 we consider a third example of the separation 
potential for the flapped airfoil which we compare with the above results. The extension of 
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the separation potential concept to 3-D wings in Section HI has the possibility of becoming 
a very useful design tool. 


Minimum section drag 


We turn now to the high Reynolds number approximation of the minimum section 
induced drag problem in Section 2. 1 . First, we introduce the asymptotic approximation of 
the drag kernel given by (2.1.35), 

Kq (jc) = -^ATo(c4x|)sinh ox 

OX 

= n sgn x for Ixl^O and <x»l (2.2.31) 

dx2^2o\x\ 6 1 1 


Use (2.2.31) in (2.1.39) to get 


1 d 
\42ltO dx 




(2.2.32) 


Now let 


F(x) = 


n 


in which case G is the solution of the universal equation 


j_d_ 
2 n dx 



(2.2.33) 


(2.2.34) 


with the minimum section drag coefficient 

= / \G{x)dx (2.2.35) 

dci 2 V 2o ' 

and the minimum drag load function 
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(2.2.36) 


l 

i(x) = c L G(x) / jG(x)dx 
-l 

and the Reynolds number independent camber shape 


l 1 

f(x) = -^\diG(i)tn\x-i\l \G(£,)di 

Z71 ' 


(2.2.37) 


We compute the solution of (2.2.34) in the following section. It is evident from (2.2.35) 
that the minimum section induced drag is 0(1 / Jo). 

Flat plate section drag 

We now compute the section induced drag in the high Reynolds number limit for 
the flat plate foil at angle of attack. The load function is 

t{x) = ^-i^- (2.2.38) 

It V 1 + X 


and the section drag coefficient of interest is 


with 




*(«!*- «0 


AT(4t|) = -^(K 0 {4t|)sinh(n) 


(2.2.39) 


(2.2.40) 


and a»l. First introduce the transformation 


x = -cos20 


£ = — cos2 (p 


(2.2.41) 


to obtain 


n/2 nil 


E = — 2 jdOcos 2 d Jdpcos 2 (p AT^ojcos 2 0-cos 2 <p j (2.2.42) 

It t\ r\ 


0 0 
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Next rotate and stretch the coordinates (0, (p) — > (s,t) with 


d = 


s-t 


9 = 


s + t 


(2.2.43) 


Also note the symmetry in s,t to get 


8 


nil s 


£ = J dsj dr(cos 2 s + cos 2 1 j AT(2crsin ssinr) 

71 0 0 


(2.2.44) 


Now use (2.2.40) and integrate by parts with respect to t. With the change of variables 


2 

T = 2<7sin s 


(2.2.45) 


we obtain finally 


£ = 1 j | Y*o(*)sinhT 

2a 

4^/2 j f sin 2 ^ — sin 2 r j sin r 

H — f ds [ dt- j A'o(2CTsin5sinr)sinh(2o'sin5sinr) 

it i i cos r- sins 


(2.2.46) 


0 0 


where we have used the definite integral. 


4 rdT 


[— ^ 0 (T)sinhr=l 

J T 


It ' T 


(2.2.47) 


For o » 1, the Bessel function can be replaced with its asymptotic expansion in (2.2.46) 
in which case we find 


and 


£ = 1 + 0 



(2.2.48) 



for (T 


(2.2.49) 


29 



This is the experimentally verified result for the section induced drag of a flat plate foil 
(Hoemer 1965, p. 7-3). We have derived this well-known Reynolds number independent 
result with a complete viscous theory. The approach to this limit is calculated numerically in 
Section 2.3. 


Profiled airfoil section drag 


The present theory was derived formally for a zero thickness foil. When the Reynolds 
number is large the load distribution becomes identical to the classical inviscid distribution 
obtained with the Kutta condition and we have shown that the section drag can be obtained by 
tilting the resultant lift by the geometric angle of attack. Another point of view is that we have 
given up all "leading edge suction" a finite force that acts at the point leading edge in the limit of 
the classical inviscid theory as the thickness tends to zero It is natural to ask whether the 
present theory can yield any useful information about the nature of the leading edge suction 
force for a profiled leading edge. To investigate this question we consider the load function 


t(x) = 


C, 

2 n \ + x + £ 


(2.2.50) 


where e is the foil nose radius referenced to the foil chord. This load exhibits the usual peak 
near the leading edge but is otherwise regular. Following the above asymptotic analysis for the 
zero thickness case we obtain the following expression for the section drag 


with 


cl AS 1 
C ° In V2 n% yfae 



^(l— T 4 ) 

(l + r^l-r 2 ) 2 ) (l+x^l + r 2 ) 2 ) 


(2.2.51) 


S = Jdx x 2 (l + 2x 2 JT(x) = 0.3366 (2.2.52) 

0 


Define the Reynolds number based on the nose radius 
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(2.2.53) 


Then 


r Uo^ = 4 U?£'!2L = 4(J£ 
v 4 vc 

cl 85 1_ 

° ln V2 ifa 

.. cl 0.342 
_ In ^Rfl 


(2.2.54) 


which indicates that the section induced drag scales with the reciprocal square root of the 
Reynolds number based on the nose radius. 

The form of (2.2.54) and the limiting form of (2.2.49) suggest the composite 
asymptotic relation 

c 2 L 1 

Q — ~ ■ 

In 1 + Q*\j-Rh 

where the constant a is determined from (2.2.54) to be 

a = 1/.342 = 3 

From data for the roughened NACA 0006, 0009 and 00012 (Abbott and von Doenhoff 
1958) we infer the experimental value 


(2.2.55) 


(2.2.56) 


“ex P = u -^ (2.2.57) 

or 1/12 the theoretical value. The two composite relations are plotted in Figure 2-5. The 
composite section drag relation obtained with the linear viscous theory is about an order of 
magnitude smaller than the experimental value for practical airfoil sections. In fact it is 
comparable in magnitude to the section drag in the so-called laminar flow bucket (see data in 
Appendix IV of Abbott and von Doenhoff). The missing element of the theory is the presence 
of turbulent flow over the section or the turbulent enstrophy 

D^=-jj-\\6>'(dV (2.2.58) 
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where cb' is the turbulent component of the vorticity. The very smooth quadratic variation of 
the section drag with the section lift suggests that the turbulent vorticity scales directly with the 
section cl so that the general form of our composite relation seems quite plausible. The 
parameter a simply undergoes a transition at some value of cl or some value of Rn in 
Figure 2-5. The theoretical curve may be considered a lower bound on the section induced 
drag if laminar flow can be maintained over most of the section. This is an area of intensive 
current research (see Holmes 1990). 



Figure 2-5. Theoretical (laminar) and experimental (turbulent) section induced 
drag for profiled airfoil sections. 
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Parabolic camber section drag 


For any load distribution that is not singular at the leading edge or trailing edge the high 
Reynolds number section drag can be computed with (2.1.31) where we can replace the kernel 
with the asymptotic approximation (2.2.31). The formula is 


1 

D A^lna 


\dxtXxy\d$ m 
-1 -1 


sgnjx-^) 

V I* -51 


(2.2.59) 


where the derivative £'(x) may be singular. To illustrate this formula we compute the section 
drag due to the elliptic load distribution 

l(x)='^Ji^7 (2.2.60) 

n 


that corresponds to parabolic camber. Break up the integral in (2.2.59) into 2 parts x < £ and 
x > ^ and interchange the order of integration to get 


Cn — 


n 2 *j2no 


\d | Vi 7 ? 


l-l 


[ dx - r — -\dx I * 

t V( 1 -*)(*-£) i ^(l-x 2 )($-x) 


(2.2.61) 


In the first integral we let 


x = 1 - 2k 2 sin 2 (p 0 <<p< n/2 



and in the second integral 


Then 


x = 1 - 2&2 sin 2 <p 0<<p<n/2 



CD 


2 cl 


i 

jdtjl-f (lEikO-Kiki)) 
-1 


(2.2.62) 


(2.2.63) 


(2.2.64) 
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Finally set 


£ = cos2<p 


(2.2.65) 


to obtain 


16 2 nl 2 

Cq = tjx ^ fd^sin 2 ^cos 2 <p[2£(sinp)-/ir(sin<p)] 

* Va j 0 

= 0.1313 c 2 / Vcr parabolic camber (2.2.66) 

In section 2.3 we show by numerical calculation that the ideal minimum section drag is 

c D . = 0.1268 c\i4o (2.2.67) 


The result for parabolic camber is only about 4 percent greater than the ideal. The load and 
camber distributions are compared in Section 2.3. 

2.3 Numerical examples 

The basic algorithm 

To calculate the solution of the load equation (2.1.15), the universal minimum drag 
load equation (2.1.39) or its high a asymptotic approximation (2.2.34), we are faced with the 
need to approximate the integral 

C(x) = \dt /($) -f F(x - ft (2.3.1) 

In i dx 

Also, C(x) is needed to evaluate the section induced drag (2.1.31). The function Fix) for 
the various cases is defined as follows: 
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F(x) = ln\x\ + e ox K 0 (a\x\) 
= AT 0 (a1x|)sinh(ix: 


load equation (2.1.15) 

universal load equation (2.1.39) or the 
section induced drag (2.1.31) 


sgn x 

Vm 


asymptotic universal load (2.2.34) or 
section drag (2.2.59) 

(2.3.2) 


To approximate the integral for all cases we consider f(x) to be constant over intervals 
x~ to x* where x n is a control point for each interval. Then 


n=\ y~ ’ 

n 

1 N 

= JZ [ F (* -X~)-F( x- *+)] 


(2.3.3) 


Next we evaluate C(x) at control points x m to get 


N 


with 


n= 1 

K mn = [ F ( X m ~ x n)~ F ( x m ~ x n )] 


(2.3.4) 


(2.3.5) 


The control points and break points are chosen according to the cosine distribution 


x = -cos 6 


(2.3.6) 


with equal spacing on 0 . We have 

d n = (n- \/2)n/ N, 
x n = — cos e n . 


=e n ±n/2N 
x^ = -cos ( 6 n ±ji/ 2 N) 


(2.3.7) 
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The discrete kernels for the various problems of interest are summarized below: 


Load equation (2.1.15) 


N 


^ Kmn L, n ~ 

n = 1 


m 


K L = — 
mn In 


In 


fa. 


x m~ x n 


m = l,...,N 


N 


c L = Jj sin °n 

’ n=l 


High a approximation ( 2.2.15 ) 


k L = — 

™ In 


In 


Xm-Xn 


X — x + 


K 

2 G 


H{x m -x n ) H[x m -x n ) 


Universal load equation (2.1.39) 


N 




n=l 


^ mn x n ^ sinh cr|x OT X n j x n jsinh (yX m - 


Also 


1 ir N 

||F|| = J F(x)dx = — ^F n sm0 n 

-1 V n=l 


so that 



( 2 . 3 . 8 ) 


( 2 . 3 . 9 ) 


( 2 . 3 . 10 ) 

4 )] 


( 2 . 3 . 10 ) 
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(2.3.11) 


c d =4/H 

i«=c L Fj\n 


a 


= £l_ 


N 

I kLf» 


n=\ 


where K is given by (2.3.8). 

High a approximation (2.2.34) 

I«A-i 

n - 1 


*mn 


1 

sgn(x m ~x n ) 

sgn(x m -x^) 

In 

u 

x m x n 

V 

i 

1 


N 


l G ll = ^Z G n sin0 n 

n = 1 


co^ 4 m 

<n = ^C„/||CI 

with /iT^, given by (2.3.9). 

Section induced drag (2.1.31) 

The inner kernel in (2.1.31) is precisely At£,. Thus 


c ° m 2 LP mSineA " 


N 


iWi = ^Z^ sin0 « 

v n = 1 


(2.3.12) 


(2.3.14) 


(2.3.15) 


(2.3.16) 
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High a approximation (l regular) 


If the load is regular at the leading and trailing edges we can use the asymptotic 
approximation of K ^ in (2.3.14) to get, 


C D 



N 

1 


m = 1 


N 

X s ^ n ^m^mn^n 


n = 1 


(2.3.17) 


with given by (2.3.13). 

All of the following numerical results were obtained with the above algorithm, except 
for the load and separation potential of the airfoil with flap. The modifications needed to 
handle the log singularity are discussed in Appendix A. 

Constant Angle of Attack 

The load distribution l(x) and the induced drag angle ccj(x) are plotted in Figures 2- 
6a and b for o = [.001, .1, 1.0, 10, 10000]. The load evolves from the low Reynolds 
number asymptotic form (see (2.2.5)) to the high Reynolds number form (see (2.2.19)) with 
the Kutta condition (cl = 2na). The induced drag angle a d(x) tends to a constant at low 
Reynolds number which indicates that the low Reynolds number section induced drag is locally 
proportional to the load. At high Reynolds number the induced drag angle becomes more and 
more concentrated near the leading edge. This reinforces our asymptotic result (see (2.2.48)) 
for the flat plate section drag, where we showed that the leading edge was the dominant 
viscous contribution. Also, see the summary results in Figure 2-9 below. 

Parabolic camber 

The load distribution and the induced drag angle for parabolic camber arc plotted in 
Figures 2-7 a and b. At low Reynolds number, the load tends to a nearly non-lifting form with 
fore and aft anti- symmetry, compared to the highly lifting angle of attack case. At high 
Reynolds number we again recover the elliptic load distribution that is obtained with inviscid 
theory and the Kutta condition. The induced drag angle is nearly linear over the section at low 
Reynolds number and virtually constant (nearly zero in Figure 2-7b) at high Reynolds number. 
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Induced drag angle 


1 




Figure 2-7a. Load distributions: Parabolic camber; o = [0.001 to 10000]. 
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Figure 2-7b. Induced drag angle: Parabolic camber; o = [0.1 to 10000] 
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We will see below that the high Reynolds number results for parabolic camber are very close to 
those for the minimum drag solution. 

Minimum drag solution 

The load distribution and the local angle of attack or negative slope of the camber shape 
are plotted in Figures 2-8a and b. The section induced drag angle at minimum drag is a 
constant for any Reynolds number as we discussed in the derivation of these solutions. This 
important result is analogous to the 3-D minimum drag result wherein the downwash or 
induced angle of attack is constant across the span. The load distribution at minimum drag and 
a — » 0 is identical to the solution for constant angle of attack. From Figure 2-8b we see that 
a = constant is the low Reynolds number camber shape in accordance with Batchelors theorem 
(see page 21). At high Reynolds number, the minimum drag solution tends to a more elliptic 
load distribution with a shape that has logarithmic cusps at the leading and trailing edge. 



Figure 2-9. Comparison of minimum drag and parabolic camber loads and 
shape for very high Reynolds number a. 
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We compare the minimum drag solution with the parabolic solution in Figure 2-9. 
These results were computed with the asymptotic formulation of the minimum drag problem 
(see (2.3.12) through (2.3.14)). The minimum drag shape has a steeper slope at the leading 
and trailing edges than the parabolic camber shape. The load distribution is also more 
concentrated at the edges. The limiting values of the section induced drag parameter 

c D / c 2 l are 0.1268 and 0.1313 respectively for the minimum drag and parabolic camber 
configurations — about a 4 percent difference. 

Lift and drag summary 

We summarize the results for all three of the above configurations in Figures 2- 10a and 
b over the Reynolds number range o = 10'^ to lCH. Each configuration amplitude (alpha or 
the camber height) is set to yield a unit section lift coefficient for a — » <». The constant alpha 
lift coefficient is singular for a — » 0 in accordance with our asymptotic result (2.2.7). For 
parabolic camber it tends to a constant (= 2.0). The ideal section lift is unity over the full 
range of Reynolds number. The section induced drag for constant alpha tends asymptotically 
to the ideal minimum drag for o 0 and to the flat plate limiting value (1/2 ti = 0.159) for 
ct — > «». For parabolic camber the section induced drag becomes infinite at low Reynolds 
number and tends asymptotically to zero as 0.1313/Va for high Reynolds number which is 
approximately 4 percent greater than the high Reynolds number minimum induced drag 
(.1268/Vcr). A curious feature of the minimum drag curve is that the section induced drag, 

c d I Ci is a maximum for a = 1.0. 

2-D Separation potentials (G °°) 

We conclude our presentation of 2-D numerical results with a comparison of the 
separation potentials for three camber configurations: 

1 . Angle of attack (9 degree) 

2. Parabolic camber (8% chord) 

3. Flap foil (20% chord at 16 degrees) 

All configurations are for a section cl = 1.0. The inviscid load distributions with Kutta 
condition are plotted in Figure 2-11 for convenient reference. 
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Section lift coefficient 



Sigma 

Figure 2-10a. Section lift coefficient versus Reynolds number o for three geometries. 



Sigma 

Hgure 2- 10b. Section induced drag versus Reynolds number c for three geometries. 
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Figure 2-11. Inviscid (Kutta) load distributions for constant angle of attack 
parabolic camber and 20% chord flap. ’ 

The normalized viscous induced camber, Va/ v , and the normalized separation 
potential, Vo 1 a v , in Figures 2- 12a and b are calculated with the high Reynolds number 
numerical approach outlined in Appendix A. The results for angle of attack and parabolic 
camber are identical to results obtained with the analytic formulas (2.2.24) and (2.2.28). 
While the viscous induced camber shape indicates a separation problem at the flap hinge line 
the best overall indicator is the separation potential in Figure 2- 12b. For case 1, constant angle 
of attack, the potential is positive over the entire chord with a logarithmically infinite value near 
the trailing edge. A very careful asymptotic analysis of the leading edge region for case 1 (not 
included herein) indicates an exponentially small region where V<r « v is singular and 

positive. If we interpret a positive value as an indication of separation on the suction side and a 
negative value as an indication for separation on the compression side, case one is most 
susceptible to separation on the suction side near the trailing edge and near the leading edge. 


45 




Normalized viscous induced camber 



Figure 2- 12a. Comparison of normalized viscous induced camber distributions. 



Figure 2- 12b. Comparison of normalized separation potentials. 
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The relatively small value at V<T « v over most of the section may indicate reattachment of the 
familiar leading edge separation bubble. 

The parabolic camber shape is susceptible to separation on the compression side of the 
leading edge and on the suction side near the trailing edge. Also, the trailing edge has a 
logarithmically singular separation potential and is probably more susceptible than the leading 
edge. For the same lift coefficient, separation on the parabolic shape would occur before that 
on the constant alpha case. 

The results for a foil section with flap are particularly interesting. The potential for 
separation near the leading edge is small but shifts to the compression side compared to the 
constant angle of attack case. Near the hinge line separation is indicated on the suction side 
downstream of the hinge line and on the compression side upstream of the hinge line. The 
trailing edge again exhibits the usual logarithmically singular potential for separation. We point 
out that these results are for a flap setting of 16 degrees which according to inviscid theory with 
Kutta would yield a unit lift coefficient. The strong indication and likely occurrence of 
separation near the hinge line is probably responsible for the reduced effectiveness of flaps. 
We suggest that the separation potential could perhaps be used as a theoretical tool for 
correlating experimental results on flap effectiveness with Reynolds number. 

In conclusion we remark that no attempt was made to compute the section drag of the 
flap foil configuration herein. It is intuitively obvious however that large induced drag angles 
will result near the hinge line which together with the logarithmically singular load will result in 
a section drag penalty. These Reynolds dependent effects can be computed with the present 
theory. 
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III. Complete Viscous Theory of the 3-D Thin Wing 


3.1 Basic Analysis 

The remarkable feature of the present theory is that the entire program can be carried 
out for a 3-D thin wing configuration. Most of the arguments follow from the 2-D case and 
will not be repeated in detail. Some features are quite different and will be noted. Some of 
the results derived herein may also be found in the AIAA lecture notes (Yates 1990). 

The " belt sander" concept is again used to remove the primary boundary layer. 
Using vector notation, the perturbation viscous problem and boundary conditions are (see 
( 2 . 1 . 2 )): 


div v = 0 


— + grad p = curl d) 

Ox 2 a 

<y = curl v diva) = 0 

On S: u = v = 0, w — — = - a(x,y ) 

ox 


(3.1.1) 


where all dimensionless quantities are defined as in Section 2. The length c is the mean 
aerodynamic chord 


= S = _d 
C ~ b~ A 


(3.1.2) 


where b is the span, S is the planform area and A is the aspect ratio. Alternatively, we 
can write the problem in terms of pressure and vorticity: 


On S: 


V 2 p = 0 


d(Q 

3x 

Q) z 


1 9 

— V d) diva) = 0 
2<t 


3v 

dx 


du 

dy 


= 0 


(3.1.3) 


which implies that a> z = 0 everywhere and we have, 
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dp _ 

1 d(Oy 

dx 

2cr dz 

dp_ 

1 d(O x 

dy 

2 a dz 


d 2 f + d t- 1 ( d(0 x (3.1.4) 

~dx^ dz 2(7^ dy dz J 


These boundary conditions are fundamental as we have discussed before (Yates 1990). 
The first two equations show how a pressure gradient along the surface produces vorticity 
transverse to the gradient — a viscous micro version of Prandtl’s wing theory. The last 
equation shows how viscosity modifies the balance of the normal pressure gradient at the 
surface. This boundary condition is responsible for the viscous correction to the usual 
inviscid load-downwash kernel as we will see below. 

With the integral momentum and energy balance discussed by Yates 1990, we 
obtain two remarkable equations for the resultant lift and drag-due-to-lift: 


Cr. — 


y 2 puis 


= jx! dydz0) *' y 


(3.1.5) 


C D 



— f dv(ci)t + a)l) 

4Ao , v x y > 


(3.1.6) 


where T is the Trefftz plane and V is the entire volume of fluid. The drag formula is 
quadratic in the two nonzero components of the interactive vorticity while the lift is linear in 
go*. This will lead directly to the well-known quadratic relation between lift and drag. 

An important difference between the 2-D and 3-D problems is that the latter has two 
components of vorticity. Fortunately, we can represent CD* and CDy in terms of a single 
function. Since 


we let 


d(O x | dtp, 
dx dy 


(3.1.7) 
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(Or 


dB 

dy 



(3.1.8) 


where B is like a vorticity stream function. With (3.1.3) we note that B must satisfy the 
diffusion equation 


fU-Lv’a 

dx 2a 


(3.1.9) 


The fundamental solutions for p and B can be used to represent the general solution in 
terms of unknown source strengths q 0 and : 


$**•*¥* 

B = jfd&V — 

s K 


with 


S-[(*-$) 2 +(y-l,) 2 + * J ] 


IK 


(3.1.10) 


(3.1.11) 


dB 

Compute the limit of p and — as z-> 0± on S to obtain 

dz 


p(x,y,0±) = T2jtq 0 (x,y) 
dB I 


dz 


= T2jzq^(x,y) 


z=0± 


Introduce the load function 


(3.1.12) 


t(x,y) = p ( 0 )-p(0 + ) = 4^7 O 


(3.1.13) 


Next we use the x momentum equation 


to get 


du ^ dp _ 1 dcpy ~ 1 d dB 

dx dx 2a dz 2 a dx dz 


( 1 dB' 

u = -\ p + — 

l 2a dz j 


(3.1.14) 
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Since u must be zero on S we obtain 


p(x,y,0±) 


J_dB 
2a dz 


lz=0± 


(3.1.15) 


which with (3. 1.12) and (3. 1.1 3) yields the source strengths q 0 and q x in terms of t 


% = 


i(x,y) 

4k 


<h=~ 


atQ^y) 

2n 


(3.1.16) 


The y momentum equation 


* = _±( p + ±f) 

dx dy\ 2a dz) 


(3.1.17) 


with (3.1.14) shows that v = 0 on S. The representations for p and B are 


oix-^-R) 




dB 


dB 


ay dx 


(3.1.18) 


Finally, we combine these representations in the vertical momentum equation, integrate 
from -oo to x and use the boundary condition on w to get the 3-D wing load equation, 


with 


JJ d$dr\ iK L (x- %,y - rj) = -a 


(3.1.19) 






d_ 

dy 


'K>- 


,o(x-R) 


) 


(3.1.20) 


The inviscid and viscous load kernels K i^(x,y) are plotted in Figure 3- la and b 
respectively on Reynolds scaled coordinates X = ax, Y = ay. Each kernel is the upwash 
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Figure 3- la. Contour plot of the inviscid load/upwash kernel function. 
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Figure 3-lb. Contour plot of the viscous load/upwash kernel function. 
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induced by a point load at the origin. Note that the inviscid kernel is much more singular 
than the viscous kernel. Also, we point out that the complete structure of the viscous 
kernel is buried on the scale of l/o inside the singular region of the inviscid kernel. 


Resultant lift 


It is instructive to derive an expression for the resultant lift in terms of t using the 
vorticity relation (3.1.5) and the integral representation of cflx- We get 


d 


^ diy n— a,21> 


But 


T(x) 

\K 


dy 


t e~ aR 2 ? 7 e °( x2+r2 ) _ 2n -a\x\ 

jdydz-j-= \dS\rdr M 

T oo 


<J e 


( 3 . 1 . 22 ) 


so that for x downstream of any point on 5 we obtain 

C L=^ tdXiy 


(3.1.23) 


which is the correct dimensionless representation of the lift coefficient in terms of the 
pressure jump or load l. 


Drag-due-to-lift 


One of the most remarkable features of the 3-D viscous wing theory is that the drag 
coefficient (3.1.6) can be evaluated explicitly in terms of the load function just as it was in 
the 2-D problem. The derivation is given in detail in Appendix B. Here we record the final 

result 


C D = JH*'! ; 1) 


ZtcA' 


dy 


dr 1 


~ th JH'lr Kd{x s- y - n) (3 - 124) 


dl 

dx 


dl 

'H 
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with 


Xo{x,y) - tn\y - Ti\ + K£{x,y) 

^£(JC.y) = |[fi(o(«-H)) + £l(o(J? + W))] (3.1.25) 

00 -t 

£ 1 (z) = Jdr^— (3.1.26) 


the standard exponential integral. This is perhaps the single most important result 
presented in this report. The drag-due-to-lift is a correlation of the spanwise and chordwise 
gradients of the load in the viscous induced drag kernel Kd — not too surprising since the 
load gradients are the source of all vorticity which contribute to the drag. This important 
induced drag kernel is plotted in Figure 3-2 with respect to Reynolds scaled coordinates 
X = gx, Y = Gy to reveal the viscous structure. 
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Figure 3-2. Contour plot of the viscous induced drag kernel function. 
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If we drop all terms that depend on the Reynolds number a we obtain one 
representation of the classical induced drag that depends only on the spanwise load 
gradients; i.e.. 


Q> 


cltxsic&l 




(3.1.27) 


with the span loading 

L= jl(x,y)dx (3.1.28) 

c(y) 

The "inviscid" logarithmic kernel in (3.1.27) cooresponds to the horizontal lines in the plot 
of Figure 3-2. For the elliptic span load distribution it is easily verified that (3.1.27) yields 
the well-known minimum drag result 


Q>, 


'daiuol ^4 


(3.1.29) 


independent of the details of the chordwise load distribution and the Reynolds number. We 
re-emphasize the point that we have derived this well-known result via the formula (3.1.6) 
that is unquestionably of viscous origin! All drag by whatever name we may give it is of 
viscous origin. 

Another immediate result can be obtained with (3.1.24) if we assume that the span 

becomes infinite and — = 0. Then the first integral is zero and the y integral of Ke can 

dy 

be evaluated to obtain our previous result (2.1.31) for the section drag due to lift. We have 
already demonstrated the power of the 2-D theory in establishing the optimal chordwise 
load distribution to minimize section drag. Also we derived a Reynolds independent drag 
component for the flat plate at angle of attack— the leading edge suction. The formula 
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(3.1.24) combines all of these features in a single expression for the drag-due-to-lift that is 
Reynolds number dependent. 

Comment on Greene’s "viscous induced drag " 

G.C. Greene (1988) proposed a modification of the traditional approach to the 
minimum induced drag problem that in principle uses the fundamental formula for drag in 
terms of enstrophy and the vorticity production formulas in terms of pressure gradients. 
Instead of minimizing the drag, however, he minimized the streamwise vorticity production 
by minimizing the square of the span load gradient which he incorrectly claims is the drag. 
The resultant span load is certainly not elliptic but the drag is also not a minimum. 
Furthermore, Greene's formulation does not reveal any of the coupled spanwise/chordwise 
structure of the viscous induced drag kernel in Figure 3-2. 

Minimum drag-due-to-lift 


The analysis of minimum drag and optimal wing loading is straightforward with the 
above formula for Cq. We construct the functional on t ^ 



(3.1.30) 


where X is a Lagrange multiplier for the constraint on C\_. The variation of (3.1.30) 
yields the Euler equation for optimal wing loading. Eliminating the Lagrange multiplier we 
obtain the following results that may be compared with those in Section 2.1: 

K D (x-^y-ri) = \ (3.1.31) 

s 

with 
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d 1 d 2 

■+ 


V 2 = ■ 

r = ^c 2 

t(x,y) = lAC L F(x,y)l\F\ 

Ill’ll = JJ ^ forty 


(3.1.33) 


The optimal local angle of attack or camber slope distribution is given by 

a(x,y)~~jjd %dl\l-K L 


(3.1.34) 


where Kl is given by (3.1.20). For a given planform the above theory will yield the 
optimal load distribution and camber shape both spanwise and chordwise. Classical 
minimum drag theory or Greene's modification can only yield information about the span 
load distribution. 


Separation potential 


By inspection of the load-downwash relation (3.1.19), we observe that the inviscid 
and viscous parts of the kernel function Kl are separable if we introduce some device for 
omitting the neighborhood of the stronger singularity of the two individual parts (e.g., 
Hadamard finite part). Thus, we can define a 3-D viscous induced angle or separation 
potential a v . We write 


j-§d$dr\ l K my =-{a+ a v ) 


(3.1.35) 


with 


a, 


(3.1.36) 


and 




(3.1.37) 
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The usual inviscid kernel is, 


Ku, 


|.I-4.[lf 1+ £Y 

dx R c?y[yv RJ_ 


(3.1.38) 


Our main interest here is in the separation potential a v . We introduce a doublet strength 
A<p such that 


/ = 


dA(p 

dx 


(3.1.39) 


where Acp is zero on the leading edge of S and takes on the value A<p e at the trailing 
edge where 


A(p e = - j £(x,y)dx = -L(y) (3.1.40) 

c(y) 

Integrate by parts in (3. 1.36) and use the above results to get 


a v (x,y) = } dr] V 2 JJ d$di] A <p 


,o(x-$-R) 


(3.1.41) 


where K‘ is the viscous kernel (3.1.37) for source points on the trailing edge 
(£ =*eCn))- Explicit examples of the separation potential are given in Section 3. For high 
Reynolds number, it provides a powerful new tool for assessing the practicality of 
planform geometries and selected camber distributions. 


3.2 High Reynolds number results 

At present the status of the 3-D theory is incomplete. Many examples remain to be 
calculated and derived asymptotically to complete the program. Nevertheless, we present 
some key results at high Reynolds number that illustrate the utility of the viscous analysis. 

When a » 1, the various viscous kernels in Section 3.1 become highly banded in 
the spanwise direction. For example, consider the viscous kernel (3.1.37) or the integrand 
in (3.1.41). The object 

e o(x-$-R) 

R 
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decays exponentially for aly-Tjl » 1. In fact the relation a(x-^-fl) = C defines a 
viscous parabolic zone of influence that opens downstream of the source point (£,r|). 
Thus the object has the behavior of a delta function and we write 

-oR 

= CS(y) (3.2.1) 

R 

with 

00 a ~°R 

C= \dy- = 2K 0 (o\x\) (3.2.2) 

* R 


which reminds us of our previous 2-D results. With the above we can write 

K v =-2£(y)JU“AT 0 (oM) 

X 

-28'\y) \d^K 0 (a\^\) (3.2.3) 


We further note that for x < 0, K w becomes exponentially small on the scale of 1/a. Thus 
we can neglect the first term in the separation potential (3.1.41) and write 

a,(.x,y) = V 2 |J did^Ave^Ko^ - 1, |)5(y - n) (3.2.4) 


Finally, we replace the Bessel function by its asymptotic approximation to get 


a v (j:,y) 


1 

2V2 KG 


V 2 - jdi 
*o(y) 




(3.2.5) 


where x Q (y) is the leading edge coordinate. When the function a v is small, the solution 
of the inviscid problem with trailing edge Kutta condition is a valid solution of the load 
problem. However, we will see that a v can become quite large near the trailing edge and 
at lateral edges that are aligned with the flow. This function can be computed routinely 
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with the output of a doublet lattice solution of the wing lift problem. It is the 2-D planform 
Laplacian of the chordwise Abel transform of the doublet strength Atp. 

Prandtl’s relations 

The differential relationship between shed stream wise vorticity and span loading is 
well known in classical wing theory. We have already noted that the first two of Equations 
(3.1.4) are micro versions of Prandtl’s early wing theory that follow directly from the 
Navier Stokes equations. With the representation of vorticity in (3.1.18) we can derive 
another intermediate load vorticity relationship that illustrates Prandtl's formulas. We 
define the local circulation components 


But 


with 


Yx = f a x dz Yy= J coydz 

-OO -OO 

^ OO QQ 

= iS Bd2 = -iS Bdz (3.2.6) 

-OO -OO 

~jBdz = -^-jfdl;dvle , X x -V 

= ~Jjd4dr,< e ‘’^%(aR) ( 3 2 . 7 ) 

s 

^ = [(* ~ 's) 2 + Cy - n) 2 ]^ (3.2.8) 


Now for ct»1, Ko(oR) has delta function properties with respect to y. Thus 

KQ(oR) = C8(y-T]) (3.2.9) 

with 

OO 

C= J> 0 (oK)dy 

— OO 

= ! e -°\ x S\ (3.2.10) 

<T 

and 
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Thus 


dz = -\\d^drit(^,ri)e 0 ^ x 

— oo S 

= -J dti(S,y) (3.2.11) 

x o 


Y x =-^] d ^(M 


*0 


= -— 
dy 

Y y = t{*>y) 
= 0 


x< x e 


x>x e 

x<x e 

x>x e 


( 3 . 2 . 12 ) 


which are the well-known wing theory relations. For x > x e (the trailing edge) the 
spanwise circulation is zero and the streamwise circulation is the negative of the span load. 

Drag-due-to-loads 

Next we consider the kernel K E (see (3.1.25)) that appears in the drag-due-to- 
loads and in the minimum drag problem. It too has delta function properties that we want 
to exploit. Consider 

^(©(*±1x0)2^500 

with 

C* = ]dy E^ct^+M)) (3.2.13) 

— oo 

Integrate by parts to get 



— oo 

= 2|x|/> jc| (/i: 1 (alx|) + AT 0 (°M)) (3.2. 14) 
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For a » 1 we obtain the asymptotic results 


and 


C + =0 



^E(x,y) = 2^j^ ^/ixj S(y) (3.2.16) 


Substitute (3.2.16) into (3.1.24) and rewrite after integration by parts as 
Co= - 

b b 

+ 4 Xjin& N ^ £ ( x ' y ) v2 \ d Z l (^y)4 x ~^\ (3-2- 17) 

S c 

where c and b denote integration over the local chord c(y) and total span b. The above 
result has been derived somewhat formally with little regard for the singular behavoir of 
t(x,y) near free edges. We have already shown with our 2-D analysis that the leading 
edge singularly leads to an order one contribution to Co — the leading edge suction or 
absence thereof. Singularities at the wing tip will also lead to a singular and perhaps order 
one contribution to Co- Further analysis will be required to verify or refute this suspision, 
but the second derivative with respect to y in (3.2.17) is a strong indication that such is 
the case. 

Minimum drag problem 

The problem for the minimum drag load distribution can also be reduced by using 
the above asymptotic approximation for Ke. We write (3.1.31) in the form 
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±£ r jj li l;dT,F($,T,)tn\ > -4 

s 

v 2 • JJ d$dri F(£, 77 ) • K e (x - $,y - tj) = 1 (3.2. 18) 


The second term is the new feature of the viscous theory. Without this term we recover the 
classical equation for the minimum induced drag span load distribution 


J_ d2 

4 it dy 1 


A 

jdi] F(ri)ln\y--n\ = l 

-A 


(3.2.19) 


with the elliptic load solution 


F(y)= $F(x,y)dx = 4^A 2 -y 2 \F\ = 2jcA 2 

c(y) 

and famous formula for induced drag: 

^2 

C D = — (classical) (3.2.20) 

itA 

If the span becomes infinite with uniform chord and we neglect all spanwise derivatives, 
we recover the 2-D minimum section drag problem that was discussed in Section II. 

A plausible method of attacking the combined problem is the following. Since the 
two integral terms in (3.2.18) are of different asymptotic order we might seek a solution 
F{x,y) such that each term individually is a constant; i.e., 

f(7]) ^ 3 '" T, l = C ° 

* b 

^ ^ = Cl (3.2.21) 

C 


With 


64 


(3.2.22) 


C I Q — 1 

c ° v? 

the global span loading F(y) will be elliptic with small 0{\ / Vtf ) Reynolds number 

corrections. Away from the tip region F(x,y) will tend to the 2-D result obtained 
previously. A more careful numerical analysis of the second of equations (3.2.21) should 
in principal yield the optimal load distribution near the tip. 

We conclude with the remark that the optimal load equation can be expressed 
compactly as the 2-D Laplacian of a single integral expression; i.e., 

V 2 JJ d$dr\ F(£, 7 i)K d (x -Z,y-ri) = \ (3.2.23) 


where Kq is the viscous induced drag kernel (3.1.25). This result follows from the fact 
that the first term in (3.2.18) is independent of x. 

3.3 Example of the separation potential 

To conclude our present treatment of the 3-D wing we derive an explicit formula for 
the high Reynolds number separation potential (3.2.5) and illustrate the result with a 
numerical example for the rectangular wing. It is convenient to work with the load form of 
(3.2.5) which we obtain with a single integration by parts. The result is 

. x 

a v (x,y) = - -j^= V 2 J d£ /(£,y) jx-% (3.3. 1) 

x 0 

where o » 1. In general we would implement this relation by solving for the inviscid 
load or doublet strength Acp with a trailing edge Kutta condition. Then calculate a v (x,y) 
numerically to obtain a map on the surface where separation can be expected — a post 
processing operation with a doublet lattice solver. Below we assume a general form of the 
load distribution and derive an analytic expression for a v . 

Consider a load distribution of the form 
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(3.3.2) 


/U.y)=-7 

Jt (x e -x 0 )\x-x 0 

where x 0 (y),x e (y) denote the leading and trailing edge of the wing respectively and L(y) 
is the span load, i.e.. 


Also 



(3.3.3) 


(3.3.4) 


where A is the aspect ratio. Substitute (3.3.2) into (3.3.1) and introduce the change of 
variables 


=xq +2(x e - x 0 )k 2 0<<p<7i/2 

k 2 = X (3.3.5) 

We obtain 

a ' ix ' y)= -^kz* 2 { L(yWk) \ (3 ' 36) 

with 

Uy) = L(y)(x t -x 0 )% 

J(k) = (l + k 2 )E{k) - (l-k 2 )K(k) (3.3.7) 


The elliptic integrals follow from results in Appendix A. The final operation is to reduce 
the derivatives that are indicated in (3.3.6). We use the known results for derivatives of E 
and K, 
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(3.3.8) 


k 

r.-±r-£ 

k(l-k 2 ) k 

with the following formula for 

dk_ 1 
dx 2 k{x e -x 0 ) 

to obtain the following: 

dJ _ Ajk) _ 3 £ 

2(* e -*<,)* 2(* e -x 0 ) 

<? 2 / = 3 1 f£z£^ 

Ih? 2(x e -x 0 ) 2{x e -x 0 )k\ k ) 

3(E - AT) 

4(x e -x 0 fk 2 

Also 

^- = k y J' = 3kk y E 

dy ’ 

fj^T = 3(* 2 + k kyy )e + 3* * 2 

= 3(2^y + k kyy^E - 3kyK (3.3.11) 


(3.3.9) 


(3.3.10) 


where k depends on y through (3.3.5). Introduce the above derivatives into (3.3.6) to 
obtain 


a y (x,y) = 


Sn-JlKa 


L(y) 


f d 2 J + d 2 J^ 
dx 1 dp 


+ 2 L'(y)^ + L"(y)J 

dy 


or 


« v = 


JMs [(* -x„)(x- x 0 ) + 8 ( 2t ? + 


(3.3.12) 
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Rectangular wing 


For the special case of a rectangular wing, k is independent of y and. 



* 2 =— , L = V2L()0 (3.3.13) 

2 


so that 


a v = ~!M- • ££ - [d + * 2 )fi - a - * 2 )jc] 

4 jt^fnc k 2 3jtyn<j L J 


(3.3.14) 


Now assume an elliptic span load 


L(y) = -C L (l-y 2 /A 2 ) 


Yl 


(3.3.15) 


where A is the aspect ratio. Then 

L " = _£%(l -y 2 /A 2 Y* (3.3.16) 

itA ' 7 

and v > 

(3.3.17) 

The first term has a logarithmic singularity along the trailing edge (fc -» 1) while the 
second term has a much stronger algebraic singularity at the cut-off rectangular tip — a result 
of the second y derivative in (3.3.1). The separation potential scales linearly with the lift 
coefficient and inversely with the square root of the Reynolds number which is typical of 
most results obtained with the linear viscous theory. The importance of a v is in its 
variation over the planform. We have plotted contours of the normalized separation 

potential 


S(x,y) = 


C L 


(3.3.18) 
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in Figure 3-3 for A = 4. Also, we have overlaid contours of the doublet strength 

X 

A(p = -\l{$,y)d$ 

-1 

= -iSp. (l - y 2 /A 2 f^/2 + sin" 1 * + (3.3.19) 

for a unit lift coefficient. These contours are the vortex field lines which all turn and 
separate at the trailing edge in accordance with the Kutta condition. The striking feature of 
the contours of S(x,y ) is the singular behavior near the tip. A strong separation zone is 
indicated that originates at the tip leading edge and spreads inboard along the trailing 
edge — a physically appealing result 



Figure 3-3. Separation potential for a rectangular wing: A = 4. 


The tip separation zone is actually indicated in the general formula (3.3.12) if we 
note that L"(y) is the second span derivative of the product of the span load L(y) times 
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the square root of the local chord (x e -x 0 ) 112 (see (3.3.7)). This formula also suggests that 
the tip separation region could be significantly reduced in size by contouring the planform 
such that the derivatives of L(y) are regular at the tip. For example, assume that L(y) is 
locally elliptic, i.e., 

L(y) = 0(\-y / A)^ 2 y->A (3.3.20) 

and that 

x e -x 0 = 0(l-y / A) v (3.3.21) 


Then 


l+i) 

L(y) = 0{l-y/ A) 2 


(3.3.22) 


and the derivatives of L are regular if v = 1. Thus, a wing tip that tapers smoothly to a 
point appears to indicate a regular behavior of the separation potential near the tip if our 
assumption of a local elliptic loading is correct. It is a curious fact that the fins of many 
marine mammals and fish and many birds have highly pointed wing tips. Also, the wing 
tips of recent transport aircraft have evolved towards this shape (e.g., the Boeing 757 and 
767). For a discussion of skewed and pointed wing tips, see Vijgen, et al. (1987). The 
small induced drag advantage that these authors report may be secondary compared to the 
reduction of the separation potential at higher angle of attack. 
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IV. Conclusions and Recommendations 


A unified viscous theory of lift and drag of thin 2-D airfoils and 3-D lifting surfaces 
has been developed in detail. The theory is "complete" in the formal sense of Birkhoff s 
reversibility paradox stated in Section II but far from complete in terms of the opportunities 
that it opens up for future research and technical applications. We summarize the main 
conclusions of our effort and recommend a number of avenues for further research and 
development 

Summary of conclusions 

1. The lift and section drag-due-to-lift of thin 2-D airfoils can be calculated uniquely with 
the viscous theory developed herein. 

2. The high Reynolds number limit of the calculated lift (cl = 27ta) is in agreement with 
inviscid theory that uses the trailing edge Kutta condition. The high Reynolds number 
section induced drag limit (c D -c^I'Ik') is in agreement with experimental results. The 
theory shows conclusively that the section induced drag is due to viscous action at the 
leading edge. 

3. The section induced drag of profiled airfoils is of order 1 / tJRn where Rn is the 
Reynolds number based on the leading edge radius. While the scaling appears to be correct 
(see Henderson and Holmes 1989), the magnitudes are about an order of magnitude 
smaller than those reported for very thin NACA four digit series foils that operate in 
turbulent flow. The calculated levels are in closer agreement with experimental values for 
the "laminar flow bucket". The role of turbulent enstrophy and drag is discussed. 

4. The theory of minimum section induced drag is a significant new contribution to the 
theory of thin airfoils to the Author's knowledge. For the same lift coefficient, the 
theoretical drag minimum is only about four percent smaller than the drag of a foil with 
parabolic camber. The slope of the camber shape, however, is much more cusped at the 
leading and trailing edges. 

5. The minimum section induced drag corresponds to a constant value of the section 
induced drag angle. The product of the load and the induced drag angle yields the local 
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drag. This result is completely analogous to a constant downwash angle for minimum drag 
in classical 3-D minimum wing drag theory. Minimum drag results are computed and 
compared with those for the constant angle of attack and parabolic camber configurations 
over 7 decades of Reynolds number from Stokes flow to 10,000. 

6. The concept of a viscous induced "separation potential" is introduced with detailed 
calculations for constant angle of attack, parabolic camber and an airfoil with trailing edge 
flap. In the high Reynolds number limit, this new function can be computed with a simple 
post processing operation on the load or doublet strength output of a standard wing theory 
code. Nonlinear boundary layer calculations are not required. The three example 
calculations indicate separation near the trailing edge and hinge line as we would expect. 
The separation potential is a valuable design tool. 

7. The viscous 3-D wing theory is formally complete even though we have not obtained as 
many numerical examples as we would like - in particular a demonstration of the optimal 
load distribution for candidate wing planforms. 

8. The load/upwash kernel function of the 3-D viscous theory is less singular than its 
inviscid counterpart. A unique solution of the load equation can be calculated without 
recourse to a Kutta condition along the trailing edge. In the high Reynolds number limit, 
the classical inviscid solution with Kutta condition will be obtained - proved by asymptotic 
analysis. 

9. One of the most (if not the most) important results of our work is a 3-D viscous induced 
drag kernel function (see Figure 3-2). The correlation of the spanwise load gradients in 
this kernel yields the classical induced drag in the high Reynolds number limit. 
Furthermore, all of the 2-D section results can be recovered for high aspect ratio wings 
when the correlation of chordwise load gradients in the kernel dominate the induced drag. 
The interaction of spanwise and chordwise gradients at the wing tip is a topic that should be 
given immediate attention with the new theory. 

10. The 3-D viscous theory of minimum induced drag yields an equation for the universal 
minimum drag load distribution both spanwise and chordwise for a given planform. The 
use of this equation to analyze and synthesize wing tip designs is highly recommended. 
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11. The concept of a viscous separation potential is extended to 3-D with an explicit high 
Reynolds number example for a rectangular wing. Separation along the wing tip starting at 
the leading edge is predicted (see Figure 3-3). Again, the 3-D separation potential can be 
calculated with the output of a panel code without resorting to boundary layer calculations. 

Recommendations 

The viscous theory and concepts developed in this report open up a new avenue of 
thinking about the fundamental viscous forces that act on a body in a moving fluid 
medium;i.e., the lift and the drag. The theory is more than a conceptual exercise in that the 
fundamental forces can be calculated from viscous first principles without recourse to 
engineering empiricisms like the Kutta condition or the leading edge suction analogy. At 
the same time, the theory does not invalidate but rather reinforces the use of these tools at a 
deeper level. The theory also provides a solid foundation for further investigation of some 
of the more difficult problems of airfoil and wing theory that involve geometric thickness, 
finite angle of attack, turbulence and compressibility. The oscillating airfoil and wing and 
motions that are designed for propulsion are amenable to the type of analysis presented 
herein. In fact the original motivation for and effort on the viscous thin airfoil problem was 
to solve the unsteady problem. The unsteady form of the resultant axial force contains the 
drag discussed herein plus an additional term that represents the work done by the pressure 
forces. This term can oppose the drag and thus makes it possible to design optimal 
propulsion configurations. There are many possible avenues of fruitful research. We 
enumerate a few of the more important ones below. 

1. Complete the program for the 3-D wing with extensive numerical examples of optimal 
minimum drag loading and separation potential for different planform shapes and camber 
distributions. Place particular emphasis on the analysis and synthesis of wing tip shapes. 

2. Include geometric thickness and mean boundary layer displacement thickness in the 
basic viscous analysis of the 2-D and 3-D thin surface. Compare calculations with the 
wealth of experimental data - in particular separation data. Quantify, if possible, the 
concept of separation potential by direct correlation with experimental data. 

3. Develop the entire viscous program for small amplitude oscillatory motion of finite 
thickness foils and wings. Drag and separation potential due to unsteady motion can be 
calculated. 
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4. Include compressibility in the viscous formulation. The viscous origin of the 
fundamental forces should become even more important as speeds become transonic. 

5. Develop the minimum drag/maximum thrust analog of the theory to aid in the design of 
optimal propulsion configurations. 

6. Develop a rational theory of drag due to turbulent enstrophy and incorporate it into the 
section induced drag theory. 


74 



References 


Abbott, Ira H. and von Doenhoff, Albert E. 1959 Theory of Wing Sections Dover 
Publications, Inc., New York. 

Abromowitz, M. and Stegun, I. A. eds. 1964 Handbook of Mathematical Functions (2nd 
Printing ) National Bureau of Standards, Washington, D. C. 

Batchelor, G. K. 1967 An Introduction to Fluid Dynamics Cambridge University Press, 
p. 227. 

Birkhoff, Garrett 1955 Hydrodynamics - a study in logic fact and similitude Dover 
Publications, Inc., New York. 

Brown, S. N. and Stewartson, K. J. 1970 Trailing Edge Stall Jour. Fluid Mech., Vol 42, 
p. 561. 

Carrier, G. F., Krook, M. and Pearson, C. E. 1966 Functions of a Complex Variable - 
Theory and Technique McGraw-Hill, Inc., New York, Ch. 8. 

Erdelyi, A. ed. 1954 Table of Integral Transforms Bateman Manuscript Project Vol. 1 
McGraw-Hill Book Co., Inc., New York. 

Greene, George C. 1988 Viscous Induced Drag AIAA Paper No. 88-2550 Presented at the 
Sixth Applied Aerodynamics Conference, Williamsburg, VA June 6-8. (Also see 
"An Entropy Method for Drag Minimization." SAE Technical Paper Series 892344 
September 1989.) 

Groebner, Wolfgang and Hofreiter, Nikolaus 1961 Unbestimmte und Bestimmte Integrate 
Springer-Verlag. 

Henderson, William P. and Holmes, Bruce J. 1989 Induced Drag - Historical Perspective 
SAE Aerotech '89 Conference, Anaheim, CA, (September 25-28). 

Hoemer, S. F. 1965 Fluid-Dynamic Drag Published by the author, Hoerner Fluid 
Dynamics, P.O. 342, Brick Town, NJ 09723. 

Holmes, Bruce J. 1990 Prospects for Drag Reduction AIAA Professional Studies Series 
1990 Drag-Prediction and Measurement Portland, Oregon (August 23-24). 

Horne, W.Clifton., Smith, Charles A. and Karamcheti, Krishnamurty 1990 Some 
Aeroacoustic and Aerodynamic Applications of the Theory of Nonequilibrium 
Thermodynamics AIAA paper 90-3989 (October 22-24). 

McLean, J. D. 1990 Viscous Drag AIAA Professional Studies Series 1990 Drag-Prediction 
and Measurement Portland, Oregon (August 23-24). 

Melnik, R. E. 1978 Wake Curvature and Trailing Edge Interaction Effects in Viscous Flow 
Over Airfoils ATAR Conference, NASA Langley Research Center, Hampton, VA 
(March 7-9). 


75 



Rooney, E. C. 1990 The Determination of Drag from Flight Tests AIAA Professional 
Studies Series 1990 Drag-Prediction and Measurement Portland, Oregon (August 
23-24). 

Shen, S. F. and Crimi, P. 1965 The Theory for an Oscillating Thin Airfoil as Derived from 
the Oseen Equations Jour. Fluid Mech., Vol. 23, part 3, pp. 585-609. 

Southwell, R. W. and Squire, H. B. 1933 A Modification of Oseen’s Approximate 
Equations for the Motion in Two Dimensions of a Viscous Incompressible Fluid 
Phil. Trans, of Roy. Soc., Series A, 232, 27. 

van der Vooren, J. and Slooff, J. W. 1990 CFD-Based Drag Prediction; State-of-the-Art, 
Theory, Prospects AIAA Professional Studies Series 1990 Drag-Prediction and 
Measurement Portland, Oregon (August 23-24). 

Vijgen, P.M.H.W., Van Dam, C.P. and Holmes, B.J. 1987 Sheared Wing-Tip 
Aerodynamics: Wind-Tunnel and Computational Investigations of Induced-Drag 
Reduction AIAA paper 87-2481 CP, 5th Applied Aerodynamics Conference, 
Monterey, CA (August 17-19). 

Yates, John E. 1990 Prediction of Drag and Lift - a Viscous/Thermodynamic Approach 
’ AIAA Professional Studies Series 1990 Drag-Prediction and Measurement 
Portland, Oregon (August 23-24). 

Yates, John E. and Donaldson, Coleman duP. 1986 A Fundamental Study of Drag and an 
Assessment of Conventional Drag-due-to-Lift Reduction Devices NASA Contractor 
Report 4004 (September). 

Yates, John E. 1986 The Development and Application of a Unified Viscous Theory of 
Drag and Lift ARAP Tech Memo 86-12 (April). 

Yates, John E. 1983 A Research Prospectus on the Viscous Theory of Lift and Related 
Problems in Aeroacoustics ARAP Tech. Memo. 83-5 (April). 

Yates, John E. 1982 A Study of Thickness and Viscous Effects in Unsteady Three- 
Dimensional Lifting Surface Theory ARAP Report No. 463 prepared for NASA 
Langley Research Center, Contract No. NAS 1-16241 (January). 

Yates, John E. 1980 Viscous Thin Airfoil Theory ARAP Report No. 413 prepared for 
Office of Naval Research, Contract No. N00014-77-C-0616 (February). (Also see 
Viscous Thin Airfoil Theory and the Kutta Condition AIAA Paper No. 78-152, 
January 1978 and Unsteady Viscous Thin Airfoil Theory AGARD Report No. 671 
January 1979). 


76 



Appendix A 

2-D Separation Potentials 


1. Constant angle of attack (see (2.2.20)) 


n (r> = 1 t f d Z c L -2ncc Z 

injnt dx J7T7 Vi-f 


Let 


so that 


£ = -l + (l + x)sin 2 0 O<0<jt/2 


VC l-f 2 )(x-£) = V2^sin20Vl->t 2 sin 2 < l> 


_ fi+x 

It 


2 


= V2- 


d0 


Va V 1 -^ sin 2 ^ 


Thus 


with 


a - w = ' -2«MW) + 4 *uem] 


*/2 


x/2 


E(k) = J d<pjl-k 2 sin 2 0, AT(A) = f d<p /J l -k 2 sin 2 0 

0 o 

With the derivative formulas. 


E'(k) = (E-K)/k, K'(k) = 


E K dk = J_ 
k(l-k 2 ) k ’ 


we obtain 


(x) = / a ■ - * 2 ) ■ - jo ■ + 4^(£ - jo] 


(A.l) 
(A. 2) 

(A.3) 

(A.4) 
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(A.7) 
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2. Parabolic camber (see(2.2.27)) 


«,=- 
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(A. 8) 
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Now integrate by parts to get 
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(A. 11) 
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3. Foil with flap 


Consider the solution of the inviscid load equation 


-jd^-^=2a(x) 

x-% 

(A. 13) 

a(x) = /3H(x-x 0 ) 

(A. 14) 


where (5 is the flap angle and x 0 is the hinge line. The solution of the load problem with 
trailing edge Kutta condition is outlined below. The numerical algorithm for evaluating 
the separation potential follows. 

Load distribution 

The general solution of (A. 13) for any function a(x) is 


Kx) = 


Wi-x 2 




For a given by (A. 14) we evaluate 
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with 
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Use the known integral (Groebner and Hofreiter, Unbestimmte Integrale (1961), p.122) 
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To satisfy the Kutta condition at x = 1, we require 
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Separation potential 


The separation potential is given by 


with 


a v 


dx 


/v = 



(A.23) 


We interpret / v (jc) as the viscous induced camber. Since the load distribution has a log 

singularity at x = jc 0 , we split the integral (A.23) into two parts x<x 0 and x > x a . Also, 
we observe that the first term in (A.22) is the load distribution for constant angle of 
attack. Thus, we can use the results (A.4) and (A.7) with c L = 2na to get. 
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For each of the three integrals (a to b) in (A.27), we set 
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The result is 
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(A. 30) 


The above formula is used in Section 2.3 to calculate the viscous induced camber for the 
2-D airfoil with flap. The separation potential is calculated by numerical differentiation 
in each interior analytic region. The numerical algorithm was also checked against the 
analytic results for the constant angle of attack and parabolic camber cases. 
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Appendix B 

3D Viscous drag-due-to-lift 


Wc derive a 3-D formula for viscous drag-due-to-lift in this appendix starting 
with (3.1.6) of the main text with the auxiliary definitions (3.1.18): 
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Substitute (B.2) into (B.l), take the x and y derivatives inside the integral over S and 
interchange the order of integration with respect to V and S. We obtain 
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with 


X(£, 77 ) = — J dxdydz • exp(( a(2x - £)) • ext> ^ £>0 (B.4) 

where R’ is given by (B.2) with £ = 77 = 0. 

We evaluate X(£, 77 ) for the space [V: x<L] where L is some large number that 
we eventually take to infinity. This step is essential since the integral is improper in the 
far wake. Write 


B-l 



X(|, ri) = - ^ J dx exp( o(2 X - £)) • I 
— ## 

with 


i,jj^ gfcg±ga 


(B.5) 


(B.6) 


Then use the double Fourier transform to obtain the Parseval relation 
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The asterisk denotes the complex conjugate in (B.7). Next we evaluate (B.8) using 
standard Fourier cosine transforms: 
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Again evaluate (B.13) using Fourier transforms. Consider 
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and with (B.5) we get 
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Integrate (B.18) by parts to obtain 
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For very large L 
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For L » crrj 2 / 4 we obtain 


£, (r(L)) = -2 ln| + + ln(L) + constant (B.28) 

But neither the constant nor ln(L) which is really ln(x) contributes to the drag coefficient 
in (B.3). Thus we can write 

X = ln|r^ + AT £ (B.29) 

with 

M6 n) = (<J(R -\Z\)) + £, (<r(/?+|£|))l (B.30) 

and 

R = (£ 2 + tj 2 ) 112 (B.31) 


Finally, we observe that ln|r^ is an irrelevant constant in the first integral of (B.3). Thus, 

with the change of notation (£, £') -* ( x , |) and ( T] , rj') -> (y, tj ) we obtain the 
final expression for the 3-D viscous induced drag coefficient: 


with 

K d ( x, y) = ln|y| + i {£, [ c{R - 1 Jt|)] + E x [ (j(R + | jc|)]) 


(B.32) 

(B.33) 


A full discussion of this important result may be found in Section 3.2 of the main text. 
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